Vascular Endothelial Growth Factor-C Treatment Enhances Cerebrospinal Fluid Outflow during Toxoplasma gondii Brain Infection but Does Not Improve Cerebral Edema.

Cerebral edema frequently develops in the setting of brain infection and can contribute to elevated intracranial pressure, a medical emergency. How excess fluid is cleared from the brain is not well understood. Previous studies have shown that interstitial fluid is transported out of the brain along perivascular channels that collect into the cerebrospinal fluid (CSF)-filled subarachnoid space. CSF is then removed from the central nervous system through venous and lymphatic routes. The current study tested the hypothesis that increasing lymphatic drainage of CSF would promote clearance of cerebral edema fluid during infection with the neurotropic parasite Toxoplasma gondii. Fluorescent microscopy and magnetic resonance imaging was used to show that C57BL/6 mice develop vasogenic edema 4 to 5 weeks after infection with T. gondii. Tracer experiments were used to evaluate how brain infection affects meningeal lymphatic function, which demonstrated a decreased rate in CSF outflow in T. gondii-infected mice. Next, mice were treated with a vascular endothelial growth factor (VEGF)-C-expressing viral vector, which induced meningeal lymphangiogenesis and improved CSF outflow in chronically infected mice. No difference in cerebral edema was observed between mice that received VEGF-C and those that rececived sham treatment. Therefore, although VEGF-C treatment can improve lymphatic outflow in mice infected with T. gondii, this effect does not lead to increased clearance of edema fluid from the brains of these mice.

The meningeal transcriptional response to traumatic brain injury and aging.

Emerging evidence suggests that the meningeal compartment plays instrumental roles in various neurological disorders, however, we still lack fundamental knowledge about meningeal biology. Here, we utilized high-throughput RNA sequencing (RNA-seq) techniques to investigate the transcriptional response of the meninges to traumatic brain injury (TBI) and aging in the sub-acute and chronic time frames. Using single-cell RNA sequencing (scRNA-seq), we first explored how mild TBI affects the cellular and transcriptional landscape in the meninges in young mice at one-week post-injury. Then, using bulk RNA-seq, we assessed the differential long-term outcomes between young and aged mice following TBI. In our scRNA-seq studies, we highlight injury-related changes in differential gene expression seen in major meningeal cell populations including macrophages, fibroblasts, and adaptive immune cells. We found that TBI leads to an upregulation of type I interferon (IFN) signature genes in macrophages and a controlled upregulation of inflammatory-related genes in the fibroblast and adaptive immune cell populations. For reasons that remain poorly understood, even mild injuries in the elderly can lead to cognitive decline and devastating neuropathology. To better understand the differential outcomes between the young and the elderly following brain injury, we performed bulk RNA-seq on young and aged meninges 1.5 months after TBI. Notably, we found that aging alone induced upregulation of meningeal genes involved in antibody production by B cells and type I IFN signaling. Following injury, the meningeal transcriptome had largely returned to its pre-injury signature in young mice. In stark contrast, aged TBI mice still exhibited upregulation of immune-related genes and downregulation of genes involved in extracellular matrix remodeling. Overall, these findings illustrate the dynamic transcriptional response of the meninges to mild head trauma in youth and aging.

Meningeal lymphatic drainage promotes T cell responses against Toxoplasma gondii but is dispensable for parasite control in the brain.

The discovery of meningeal lymphatic vessels that drain the CNS has prompted new insights into how immune responses develop in the brain. In this study, we examined how T cell responses against CNS-derived antigen develop in the context of infection. We found that meningeal lymphatic drainage promotes CD4+ and CD8+ T cell responses against the neurotropic parasite Toxoplasma gondii in mice, and we observed changes in the dendritic cell compartment of the dural meninges that may support this process. Indeed, we found that mice chronically, but not acutely, infected with T. gondii exhibited a significant expansion and activation of type 1 and type 2 conventional dendritic cells (cDC) in the dural meninges. cDC1s and cDC2s were both capable of sampling cerebrospinal fluid (CSF)-derived protein and were found to harbor processed CSF-derived protein in the draining deep cervical lymph nodes. Disrupting meningeal lymphatic drainage via ligation surgery led to a reduction in CD103+ cDC1 and cDC2 number in the deep cervical lymph nodes and caused an impairment in cDC1 and cDC2 maturation. Concomitantly, lymphatic vessel ligation impaired CD4+ and CD8+ T cell activation, proliferation, and IFN-γ production at this site. Surprisingly, however, parasite-specific T cell responses in the brain remained intact following ligation, which may be due to concurrent activation of T cells at non-CNS-draining sites during chronic infection. Collectively, our work reveals that CNS lymphatic drainage supports the development of peripheral T cell responses against T. gondii but remains dispensable for immune protection of the brain.

Microglial STAT1-sufficiency is required for resistance to toxoplasmic encephalitis.

Toxoplasma gondii is a ubiquitous intracellular protozoan parasite that establishes a life-long chronic infection largely restricted to the central nervous system (CNS). Constant immune pressure, notably IFN-γ-STAT1 signaling, is required for preventing fatal pathology during T. gondii infection. Here, we report that abrogation of STAT1 signaling in microglia, the resident immune cells of the CNS, is sufficient to induce a loss of parasite control in the CNS and susceptibility to toxoplasmic encephalitis during the early stages of chronic infection. Using a microglia-specific genetic labeling and targeting system that discriminates microglia from blood-derived myeloid cells that infiltrate the brain during infection, we find that, contrary to previous in vitro reports, microglia do not express inducible nitric-oxide synthase (iNOS) during T. gondii infection in vivo. Instead, transcriptomic analyses of microglia reveal that STAT1 regulates both (i) a transcriptional shift from homeostatic to "disease-associated microglia" (DAM) phenotype conserved across several neuroinflammatory models, including T. gondii infection, and (ii) the expression of anti-parasitic cytosolic molecules that are required for eliminating T. gondii in a cell-intrinsic manner. Further, genetic deletion of Stat1 from microglia during T. gondii challenge leads to fatal pathology despite largely equivalent or enhanced immune effector functions displayed by brain-infiltrating immune populations. Finally, we show that microglial STAT1-deficiency results in the overrepresentation of the highly replicative, lytic tachyzoite form of T. gondii, relative to its quiescent, semi-dormant bradyzoite form typical of chronic CNS infection. Our data suggest an overall protective role of CNS-resident microglia against T. gondii infection, illuminating (i) general mechanisms of CNS-specific immunity to infection (ii) and a clear role for IFN-STAT1 signaling in regulating a microglial activation phenotype observed across diverse neuroinflammatory disease states.

Investigating the genetic architecture of eye colour in a Canadian cohort.

Eye color is highly variable in populations with European ancestry, ranging from low to high quantities of melanin in the iris. Polymorphisms in the HERC2/OCA2 locus have the largest effect on eye color in these populations, although other genomic regions also influence eye color. We performed genome-wide association studies of eye color in a Canadian cohort of European ancestry (N = 5,641) and investigated candidate causal variants. We uncovered several candidate causal signals in the HERC2/OCA2 region, whereas other loci likely harbor a single causal signal. We observed colocalization of eye color signals with the expression or methylation profiles of cultured primary melanocytes. Genetic correlations of eye and hair color suggest high genome-wide pleiotropy, but locus-level differences in the genetic architecture of both traits. Overall, we provide a better picture of the polymorphisms underpinning eye color variation, which may be a consequence of specific molecular processes in the iris melanocytes.

A large Canadian cohort provides insights into the genetic architecture of human hair colour.

Hair colour is a polygenic phenotype that results from differences in the amount and ratio of melanins located in the hair bulb. Genome-wide association studies (GWAS) have identified many loci involved in the pigmentation pathway affecting hair colour. However, most of the associated loci overlap non-protein coding regions and many of the molecular mechanisms underlying pigmentation variation are still not understood. Here, we conduct GWAS meta-analyses of hair colour in a Canadian cohort of 12,741 individuals of European ancestry. By performing fine-mapping analyses we identify candidate causal variants in pigmentation loci associated with blonde, red and brown hair colour. Additionally, we observe colocalization of several GWAS hits with expression and methylation quantitative trait loci (QTLs) of cultured melanocytes. Finally, transcriptome-wide association studies (TWAS) further nominate the expression of EDNRB and CDK10 as significantly associated with hair colour. Our results provide insights on the mechanisms regulating pigmentation biology in humans.

A UVB-responsive common variant at chromosome band 7p21.1 confers tanning response and melanoma risk via regulation of the aryl hydrocarbon receptor, AHR.

Genome-wide association studies (GWASs) have identified a melanoma-associated locus on chromosome band 7p21.1 with rs117132860 as the lead SNP and a secondary independent signal marked by rs73069846. rs117132860 is also associated with tanning ability and cutaneous squamous cell carcinoma (cSCC). Because ultraviolet radiation (UVR) is a key environmental exposure for all three traits, we investigated the mechanisms by which this locus contributes to melanoma risk, focusing on cellular response to UVR. Fine-mapping of melanoma GWASs identified four independent sets of candidate causal variants. A GWAS region-focused Capture-C study of primary melanocytes identified physical interactions between two causal sets and the promoter of the aryl hydrocarbon receptor (AHR). Subsequent chromatin state annotation, eQTL, and luciferase assays identified rs117132860 as a functional variant and reinforced AHR as a likely causal gene. Because AHR plays critical roles in cellular response to dioxin and UVR, we explored links between this SNP and AHR expression after both 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) and ultraviolet B (UVB) exposure. Allele-specific AHR binding to rs117132860-G was enhanced following both, consistent with predicted weakened AHR binding to the risk/poor-tanning rs117132860-A allele, and allele-preferential AHR expression driven from the protective rs117132860-G allele was observed following UVB exposure. Small deletions surrounding rs117132860 introduced via CRISPR abrogates AHR binding, reduces melanocyte cell growth, and prolongs growth arrest following UVB exposure. These data suggest AHR is a melanoma susceptibility gene at the 7p21.1 risk locus and rs117132860 is a functional variant within a UVB-responsive element, leading to allelic AHR expression and altering melanocyte growth phenotypes upon exposure.

Cell-type-specific meQTLs extend melanoma GWAS annotation beyond eQTLs and inform melanocyte gene-regulatory mechanisms.

Although expression quantitative trait loci (eQTLs) have been powerful in identifying susceptibility genes from genome-wide association study (GWAS) findings, most trait-associated loci are not explained by eQTLs alone. Alternative QTLs, including DNA methylation QTLs (meQTLs), are emerging, but cell-type-specific meQTLs using cells of disease origin have been lacking. Here, we established an meQTL dataset by using primary melanocytes from 106 individuals and identified 1,497,502 significant cis-meQTLs. Multi-QTL colocalization with meQTLs, eQTLs, and mRNA splice-junction QTLs from the same individuals together with imputed methylome-wide and transcriptome-wide association studies identified candidate susceptibility genes at 63% of melanoma GWAS loci. Among the three molecular QTLs, meQTLs were the single largest contributor. To compare melanocyte meQTLs with those from malignant melanomas, we performed meQTL analysis on skin cutaneous melanomas from The Cancer Genome Atlas (n = 444). A substantial proportion of meQTL probes (45.9%) in primary melanocytes is preserved in melanomas, while a smaller fraction of eQTL genes is preserved (12.7%). Integration of melanocyte multi-QTLs and melanoma meQTLs identified candidate susceptibility genes at 72% of melanoma GWAS loci. Beyond GWAS annotation, meQTL-eQTL colocalization in melanocytes suggested that 841 unique genes potentially share a causal variant with a nearby methylation probe in melanocytes. Finally, melanocyte trans-meQTLs identified a hotspot for rs12203592, a cis-eQTL of a transcription factor, IRF4, with 131 candidate target CpGs. Motif enrichment and IRF4 ChIP-seq analysis demonstrated that these target CpGs are enriched in IRF4 binding sites, suggesting an IRF4-mediated regulatory network. Our study highlights the utility of cell-type-specific meQTLs.

Astrocytes promote a protective immune response to brain Toxoplasma gondii infection via IL-33-ST2 signaling.

It is of great interest to understand how invading pathogens are sensed within the brain, a tissue with unique challenges to mounting an immune response. The eukaryotic parasite Toxoplasma gondii colonizes the brain of its hosts, and initiates robust immune cell recruitment, but little is known about pattern recognition of T. gondii within brain tissue. The host damage signal IL-33 is one protein that has been implicated in control of chronic T. gondii infection, but, like many other pattern recognition pathways, IL-33 can signal peripherally, and the specific impact of IL-33 signaling within the brain is unclear. Here, we show that IL-33 is expressed by oligodendrocytes and astrocytes during T. gondii infection, is released locally into the cerebrospinal fluid of T. gondii-infected animals, and is required for control of infection. IL-33 signaling promotes chemokine expression within brain tissue and is required for the recruitment and/or maintenance of blood-derived anti-parasitic immune cells, including proliferating, IFN-γ-expressing T cells and iNOS-expressing monocytes. Importantly, we find that the beneficial effects of IL-33 during chronic infection are not a result of signaling on infiltrating immune cells, but rather on radio-resistant responders, and specifically, astrocytes. Mice with IL-33 receptor-deficient astrocytes fail to mount an adequate adaptive immune response in the CNS to control parasite burden-demonstrating, genetically, that astrocytes can directly respond to IL-33 in vivo. Together, these results indicate a brain-specific mechanism by which IL-33 is released locally, and sensed locally, to engage the peripheral immune system in controlling a pathogen.

Meningeal lymphatic dysfunction exacerbates traumatic brain injury pathogenesis.

Traumatic brain injury (TBI) is a leading global cause of death and disability. Here we demonstrate in an experimental mouse model of TBI that mild forms of brain trauma cause severe deficits in meningeal lymphatic drainage that begin within hours and last out to at least one month post-injury. To investigate a mechanism underlying impaired lymphatic function in TBI, we examined how increased intracranial pressure (ICP) influences the meningeal lymphatics. We demonstrate that increased ICP can contribute to meningeal lymphatic dysfunction. Moreover, we show that pre-existing lymphatic dysfunction before TBI leads to increased neuroinflammation and negative cognitive outcomes. Finally, we report that rejuvenation of meningeal lymphatic drainage function in aged mice can ameliorate TBI-induced gliosis. These findings provide insights into both the causes and consequences of meningeal lymphatic dysfunction in TBI and suggest that therapeutics targeting the meningeal lymphatic system may offer strategies to treat TBI.

Massively parallel reporter assays of melanoma risk variants identify MX2 as a gene promoting melanoma.

Genome-wide association studies (GWAS) have identified ~20 melanoma susceptibility loci, most of which are not functionally characterized. Here we report an approach integrating massively-parallel reporter assays (MPRA) with cell-type-specific epigenome and expression quantitative trait loci (eQTL) to identify susceptibility genes/variants from multiple GWAS loci. From 832 high-LD variants, we identify 39 candidate functional variants from 14 loci displaying allelic transcriptional activity, a subset of which corroborates four colocalizing melanocyte cis-eQTL genes. Among these, we further characterize the locus encompassing the HIV-1 restriction gene, MX2 (Chr21q22.3), and validate a functional intronic variant, rs398206. rs398206 mediates the binding of the transcription factor, YY1, to increase MX2 levels, consistent with the cis-eQTL of MX2 in primary human melanocytes. Melanocyte-specific expression of human MX2 in a zebrafish model demonstrates accelerated melanoma formation in a BRAFV600E background. Our integrative approach streamlines GWAS follow-up studies and highlights a pleiotropic function of MX2 in melanoma susceptibility.

Meta-analysis of GWA studies provides new insights on the genetic architecture of skin pigmentation in recently admixed populations.

BACKGROUND: Association studies in recently admixed populations are extremely useful to identify the genetic architecture of pigmentation, due to their high genotypic and phenotypic variation. However, to date only four Genome-Wide Association Studies (GWAS) have been carried out in these populations. RESULTS: We present a GWAS of skin pigmentation in an admixed sample from Cuba (N = 762). Additionally, we conducted a meta-analysis including the Cuban sample, and admixed samples from Cape Verde, Puerto Rico and African-Americans from San Francisco. This meta-analysis is one of the largest efforts so far to characterize the genetic basis of skin pigmentation in admixed populations (N = 2,104). We identified five genome-wide significant regions in the meta-analysis, and explored if the markers observed in these regions are associated with the expression of relevant pigmentary genes in human melanocyte cultures. In three of the regions identified in the meta-analysis (SLC24A5, SLC45A2, and GRM5/TYR), the association seems to be driven by non-synonymous variants (rs1426654, rs16891982, and rs1042602, respectively). The rs16891982 polymorphism is strongly associated with the expression of the SLC45A2 gene. In the GRM5/TYR region, in addition to the rs1042602 non-synonymous SNP located on the TYR gene, variants located in the nearby GRM5 gene have an independent effect on pigmentation, possibly through regulation of gene expression of the TYR gene. We also replicated an association recently described near the MFSD12 gene on chromosome 19 (lead variant rs112332856). Additionally, our analyses support the presence of multiple signals in the OCA2/HERC2/APBA2 region on chromosome 15. A clear causal candidate is the HERC2 intronic variant rs12913832, which has a profound influence on OCA2 expression. This variant has pleiotropic effects on eye, hair, and skin pigmentation. However, conditional and haplotype-based analyses indicate the presence of other variants with independent effects on melanin levels in OCA2 and APBA2. Finally, a follow-up of genome-wide signals identified in a recent GWAS for tanning response indicates that there is a substantial overlap in the genetic factors influencing skin pigmentation and tanning response. CONCLUSIONS: Our meta-analysis of skin pigmentation GWAS in recently admixed populations provides new insights about the genetic architecture of this complex trait.

Cell-type-specific eQTL of primary melanocytes facilitates identification of melanoma susceptibility genes.

Most expression quantitative trait locus (eQTL) studies to date have been performed in heterogeneous tissues as opposed to specific cell types. To better understand the cell-type-specific regulatory landscape of human melanocytes, which give rise to melanoma but account for <5% of typical human skin biopsies, we performed an eQTL analysis in primary melanocyte cultures from 106 newborn males. We identified 597,335 cis-eQTL SNPs prior to linkage disequilibrium (LD) pruning and 4997 eGenes (FDR < 0.05). Melanocyte eQTLs differed considerably from those identified in the 44 GTEx tissue types, including skin. Over a third of melanocyte eGenes, including key genes in melanin synthesis pathways, were unique to melanocytes compared to those of GTEx skin tissues or TCGA melanomas. The melanocyte data set also identified trans-eQTLs, including those connecting a pigmentation-associated functional SNP with four genes, likely through cis-regulation of IRF4 Melanocyte eQTLs are enriched in cis-regulatory signatures found in melanocytes as well as in melanoma-associated variants identified through genome-wide association studies. Melanocyte eQTLs also colocalized with melanoma GWAS variants in five known loci. Finally, a transcriptome-wide association study using melanocyte eQTLs uncovered four novel susceptibility loci, where imputed expression levels of five genes (ZFP90, HEBP1, MSC, CBWD1, and RP11-383H13.1) were associated with melanoma at genome-wide significant P-values. Our data highlight the utility of lineage-specific eQTL resources for annotating GWAS findings, and present a robust database for genomic research of melanoma risk and melanocyte biology.

Loci associated with skin pigmentation identified in African populations.

Despite the wide range of skin pigmentation in humans, little is known about its genetic basis in global populations. Examining ethnically diverse African genomes, we identify variants in or near SLC24A5, MFSD12, DDB1, TMEM138, OCA2, and HERC2 that are significantly associated with skin pigmentation. Genetic evidence indicates that the light pigmentation variant at SLC24A5 was introduced into East Africa by gene flow from non-Africans. At all other loci, variants associated with dark pigmentation in Africans are identical by descent in South Asian and Australo-Melanesian populations. Functional analyses indicate that MFSD12 encodes a lysosomal protein that affects melanogenesis in zebrafish and mice, and that mutations in melanocyte-specific regulatory regions near DDB1/TMEM138 correlate with expression of ultraviolet response genes under selection in Eurasians.

A common intronic variant of PARP1 confers melanoma risk and mediates melanocyte growth via regulation of MITF.

Previous genome-wide association studies have identified a melanoma-associated locus at 1q42.1 that encompasses a ∼100-kb region spanning the PARP1 gene. Expression quantitative trait locus (eQTL) analysis in multiple cell types of the melanocytic lineage consistently demonstrated that the 1q42.1 melanoma risk allele (rs3219090[G]) is correlated with higher PARP1 levels. In silico fine-mapping and functional validation identified a common intronic indel, rs144361550 (-/GGGCCC; r2 = 0.947 with rs3219090), as displaying allele-specific transcriptional activity. A proteomic screen identified RECQL as binding to rs144361550 in an allele-preferential manner. In human primary melanocytes, PARP1 promoted cell proliferation and rescued BRAFV600E-induced senescence phenotypes in a PARylation-independent manner. PARP1 also transformed TERT-immortalized melanocytes expressing BRAFV600E. PARP1-mediated senescence rescue was accompanied by transcriptional activation of the melanocyte-lineage survival oncogene MITF, highlighting a new role for PARP1 in melanomagenesis.